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A  REFINED  NONLINEAR  ANALYSIS  OF 
LAMINATED  COMPOSITE  PLATES  ANO  SHELLS 

J.  N.  Reddy 

Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg.  VA  24061 

SUMMARY 

This  final  technical  report  summarizes  In  a  compact  form  the 
results  of  a  two-year  research  program  on  the  development  of  refined 
shear  deformation  theories  of  plates  and  shells,  and  their  analytical 
solutions.  The  detailed  results  are  reported  In  various  reports  and 
technical  papers  during  the  course  of  the  project.  A  third-order, 
nonlinear  shear  deformation  shell  theory  and  finite  element  model  that 
accounts  for  parabolic  distribution  of  transverse  shear  stresses  through 
thickness  and  the  von  Karman  nonlinear  strains  was  developed  during  an 
Investigation  sponsored  by  NASA  Langley  Research  Center  during  the  first 
year  of  this  research.  The  theory  Is  further  refined  and  extended  and 
analytical  solutions  are  developed  during  this  research.  The  most 
significant  contributions  of  this  research  are: 

1.  The  development  of  analytical  solutions  of  the  linear  theory 
for  various  boundary  conditions.  The  space-time  concept  Is 
used  to  develop  such  solutions. 

2.  The  development  of  moderate  rotation  theory  of  laminated 
composite  plates. 

3.  The  development  of  finite  element  models  for  first-ply  and 
post-first-ply  failure  analysis. 

Future  research  related  to  refined  theories  should  consider  the 
effects  of  moderate  and  large  rotations  and  the  development  of  mlcro- 
and  macromechanical  constitutive  models  for  a  better  understanding  of 
the  structural  behavior  and  prediction  of  progressive  damage  and  failure 
In  composite  structures. 

Ill 


1.  INTRODUCTION 

The  advent  of  new  composite  materials  and  their  Increasing  use  In 
various  fields  of  advanced  technology  has  generated  a  new  interest  In 
the  development  and  solution  of  consistent  refined  theories  of 
anisotropic  composite  plates  and  shells.  This  Interest  Is  due  to  the 
fact  that  the  classical  plate  theory.  In  terms  of  Its  basic  assumptions 
(l.e.  the  Klrchhoff  hypothesis),  comes  In  conflict  with  real  behavior  of 
these  new  materials.  For  example.  In  contrast  to  the  basic  assumption 
of  Infinite  rigidity  In  transverse  shear  In  the  classical  plate  theory, 
the  new  composite  materials  exhibit  a  finite  rigidity  In  transverse 
shear.  This  property  requires  the  Incorporation  of  transverse  shear 
deformation  effects. 

In  addition  to  other  shortcomings,  the  classical  plate  theory 
Involves  a  contradiction  between  the  number  of  boundary  conditions 
physically  required  to  be  fulfilled  on  a  free  boundary  and  the  number 
available  in  theory,  which  Is  to  be  consistent  with  the  order  of  the 
associated  governing  equations  (see  Stoker  (II).  The  non-fulfillment  of 
boundary  conditions  on  the  bounding  surfaces  constitutes  another  feature 
of  the  classical  theory.  In  recent  years  attempts  were  made  to  refine 
the  classical  theory  by:  (1)  Incorporating  transverse  shear  effects, 

(11)  removing  the  contradiction  which  concerns  the  number  of  boundary 
conditions  to  be  prescribed  at  each  edge,  and  (111)  fulfilling  the 
boundary  conditions  on  the  bounding  surfaces  and.  In  the  case  of 
laminated  composite  plates  and  shells,  of  the  continuity  conditions  at 
the  Interfaces  between  the  contiguous  layers.  In  addition,  the  refined 
transverse  shear  deformation  theories  can  be  used  to  model  such 
anisotropic  plates  and  shells  whose  material  exhibits  high  degree  of 
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anisotropy,  and  are  not  restricted  to  the  thinness  requirement  Implied 
by  the  classical  laminate  theory.  Another  feature  of  refined  laminate 
theories  concerns  the  adequate  Incorporation  of  the  dynamical  effects 
allowing  the  evaluation  of  the  lowest  and  higher  natural  frequencies. 

The  shear  deformation  theories  known  In  the  literature  can  be 
grouped  Into  two  classes:  (1)  stress-based  theories,  and  (2) 
displacement-based  theories.  The  first  stress-based  transverse  shear 
deformable  plate  theory  Is  due  to  Relssner  [2-4].  The  distribution 
across  the  thickness  of  the  transverse  normal  and  shear  stresses  Is 
determined  through  Integration  over  the  thickness  of  the  equilibrium 
equations  of  the  3-0  elasticity  theory.  The  associated  field  equations 
and  boundary  conditions  expressed  in  terms  of  2-0  quantities  can  be 
determined  by  using  the  variational  principles  of  the  3-0  elasticity 
theory,  or  by  considering  the  moments  of  nth  order  of  the  basic 
equations  of  3-0  elasticity  theory.  Both  methods  allow  the  reduction  of 
the  3-D  problems  to  a  2-0  equivalent  one. 

The  pioneering  work  of  the  displacement-based  theories  Is  due  to 
Basset  [5].  Based  on  Basset's  representation  of  displacement  field, 
Hildebrand,  Relssner  and  Thomas  [6]  developed  a  varlatlonally  consistent 
first  order  theory  for  shells.  The  field  equations  were  derived  using 
the  principle  of  minimum  total  potential  energy.  This  results  In  five 
equilibrium  equations  expressed  In  terms  of  five  displacement 
quantities. 

By  using  the  displacement  representation  of  Basset,  Hindi  In  [7] 
extended  Hencky's  theory  (8)  of  Isotropic  plates  to  the  dynamic  case. 
Historical  evidence  (from  the  review  of  the  literature)  points  out  that 
the  basic  Idea  of  the  displacement-based  first-order  shear  deformation 
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theory  came  from  Basset  (51,  Hildebrand,  Relssner  and  Thomas  [6]  and 
Hencky  [8] .  Mlndlln  should  be  credited  with  the  extension  of  the  theory 
to  the  dynamic  case.  The  shear  deformation  theory  of  Hencky-MIndlln  Is 
referred  as  the  first-order  transverse  shear  deformation  theory  (see 
Reddy  (9j). 

Following  these  works,  many  extensions  and  applications  of  the  two 
classes  of  theories  were  reported  In  the  literature  (see  e.g. 
Ambartsumian  [10-12],  Boal  and  Relssner  [13],  Bolle  [14],  Cheng  [15], 
Gol'denvelzer  [16-18],  Green  [19-21],  Kromm  [22,23],  Levinson  [24], 
Llbrescu  [25-30],  Lo,  Christensen  and  Wu  [31],  Medwadowskl  [32],  Murthy 
[33],  Nelson  and  Lorch  [34],  Pagano  [35,36],  Reddy  [37-40],  Rehfleld  and 
Vallsetty  [41],  Relssner  [42-47],  Schmidt  [48],  Shlrakawa  [49],  Volterra 
[50],  Whitney  and  colleqgues  [51-54],  Wilson  and  Boresl  [55],  Yang, 
Norris  and  Stavsky  [56]). 

Extension  of  the  displacement-based  theory  to  the  moderately  large 
deflections  case  is  due  to  Medwadowskl  [32]  and  the  extension  to 
laminated  plates  Is  due  to  Whitney  [51]  and  Whitney  and  Pagano  [52]. 
Schmidt  [48]  presented  an  extension  of  Kromm' s  theory  by  accounting  for 
moderately  large  deflections  (l.e..  In  the  Von-Karman  sense). 

The  second-  and  higher-order  displacement-based  shear  deformation 
theories  have  been  Investigated  by  Nelson  and  Lorch  [34],  Llbrescu  [25- 
28],  Lo,  Christensen  and  Wu  [31],  Levinson  [24],  Murthy  [33]  and  Reddy 
[37-39].  Levinson  [24]  and  Murthy  [33]  presented  a  third-order  theory 
In  which  transverse  normals  are  assumed  to  be  Inextenslblle.  The  nine 
displacement  functions  were  reduced  to  five  by  requiring  that  the 
transverse  shear  stresses  vanish  on  the  bounding  planes  of  the  plate. 
However,  both  authors  (and  also  Schmidt  [48])  used  the  equilibrium 
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equations  of  the  first-order  theory  in  their  analysis  and  so  they  are 
variational ly  Inconsistent.  As  a  consequence,  the  higher-order  terms  of 
the  displacement  field  are  accounted  for  only  in  the  calculation  of  the 
strains  but  not  In  the  governing  differential  equations  or  In  the 
boundary  conditions.  Recently,  during  the  course  of  the  present 
research,  Reddy  [37-391  corrected  these  theories  by  deriving  the 
governing  differential  equations  by  means  of  the  virtual  work 
principle.  The  theory  presented  In  [37]  accounts  for  the  von  Karman 
strains  but  Is  limited  to  orthotropic  plates,  while  that  In  [38]  deals 
with  the  small-deflection  theory  of  laminated  plates. 

Besides  the  works  cited  above,  there  exists  numerous  other  works 
pertaining  to  refined  theories  of  elastic  plates.  Thus,  In  addition  to 
the  works  of  Basset,  Relssner,  Hencky  and  Mlndlln  that  have  Inspired  the 
elaboration  of  a  great  number  of  works,  we  mention  the  symbolic  method 
advanced  by  Lur'e  [57]  which  has  generated  a  series  of  works  devoted  to 
the  refined  theory  of  the  Isotropic  flat  plates  (see  e.g.,  Westbrook 
[58,59]).  Starting  from  the  three-dimensional  field  equations  of 
elasticity,  methods  based  on  the  expansion  of  all  the  functions 
characterizing  the  plane  state  of  stress  Into  power  series  over  the 
plate  thickness  (see  Schlpper  [60]),  into  Legendre's  polynomials  (see 
Cicala  [61]),  Into  biharmonic  polynomials  (see  Teodorescu  [62]),  or  that 
resulting  by  considering  the  nth  order  moments  of  the  equations  of 
three-dimensional  elasticity  (see  Tlffen  [63],  Tlffen  and  Lowe  [64])  are 
developed. 

Another  method  which  leads  to  the  static  and  dynamic  theory  of 
Isotropic  and  anisotropic  plates.  Is  the  method  of  asymptotic 
Integration  of  the  three-dimensional  equations  of  elasticity  (see 


Gol'denvelzer  [17,181,  Relssner  [42),  and  Green  and  Naghdl  [201).  Other 
methods  are  of  mixed  character  which  use  first,  the  method  of  expanding 
the  functions  (defining  the  stress  and  strain  state)  Into  Legendre's 
polynomials,  and  then  the  asymptotic  Integration  method  (see,  for 
example,  Ponlatovskll  [65]). 

The  analysis  of  laminated  composite  plates  for  bending  and  natural 
vibration  has  received  widespread  attention  In  recent  years.  3-D 
elasticity  solutions  for  the  bending  (see  Pagano  [36]),  vibration  and 
buckling  of  simply  supported  thick  orthotropic  rectangular  plates  and 
laminates  were  obtained  by  Srlnlvas  and  his  co-workers  [66-68].  The 
Navler  solution  of  simply  supported  rectangular  plates  was  developed  by 
Whitney  and  Lelssa  [52]  for  classical  laminate  theory,  Pagano  [35], 
Whitney  [51],  Bert  and  Chen  [69]  and  Reddy  and  Chao  [39,70]  for  the 
first-order  transverse  shear  deformation  theory,  and  by  Reddy  and  his 
co-workers  [37-39,71]  for  a  refined  shear  deformation  theory. 

Papers  dealing  with  the  response  of  plates  excited  by  dynamic  loads 
of  known  time  history  have  been  less  In  number.  Sun  and  Whitney  [72,73] 
have  analyzed  the  response  of  anisotropic  plates  In  cylindrical  bending 
using  Hindi  In's  [7]  theory  to  account  for  normal  shear  stiffness.  Sun 
and  Chattopadhyay  [74]  used  the  plate  equations  developed  by  Whitney  and 
Pagano  [53]  to  analyze  a  specially  orthotropic  plate  subject  to  a  center 
Impact.  Oobyns  [75]  presented  an  analysis  of  simply-supported 
orthotropic  plates  subjected  to  static  and  dynamic  loading  conditions. 
The  vibration  frequencies  and  mode  shapes  were  then  determined  and 
solutions  for  plate  deflections,  bending  strains,  and  normal  shear 
forces  due  to  several  transient  loads  were  obtained.  Reddy  [76,77] 
obtained  the  exact  form  of  the  spatial  variation  of  the  solution  to 


forced  motions  of  rectangular  composite  plates  for  two  different 
lamination  schemes,  under  appropriate  simply  supported  boundary 
conditions.  He  reduced  the  problem  to  the  solution  of  a  system  of 
ordinary  differential  equations  In  time,  which  are  then  Integrated 
numerically  using  Newmark's  direct  integration  method. 

In  geometrically  nonlinear  theories  of  elastic  anisotropic  plates 
one  often  assumes  that  the  strains  and  rotations  about  the  normal  to  the 
midplane  are  Infinitesimal,  and  retains  the  products  and  squares  of  the 
derivatives  of  the  transverse  deflection  in  the  strain-displacement 
equations  (the  von  Karman  assumption;  see  [[32,37,78-80]).  The  full 
geometric  nonlinearity  (Implied  by  the  strain-displacement  equations  of 
nonlinear  elasticity)  In  shell  theories  was  considered  by  Naghdl  [81], 
Llbrescu  [29],  Yokas  and  Matsunaga  [82],  Hablp  [83],  and  Pletraszkiewlcz 
[84],  among  others.  Consideration  of  full  geometric  nonlinearity  not 
only  results  in  complex  equations,  but  not  warranted  In  most  practical 
problems.  On  the  other  hand,  the  von  Karman  nonlinear  theory  does  not 
account  for  moderate  rotation  terms  that  could  be  of  significance  in  the 
analysis  (especially  In  stability  problems)  of  plates  while  accounting 
for  the  transverse  normal  and  shear  strains.  The  small  strain  and 
moderate  rotation  concept  was  used  In  the  classical  theory  of  plates  and 
shells  by  Sanders  [85],  Kolter  [86],  Reissner  [87]  and  Pletraszkiewlcz 
[88],  and  In  first-order  plate  and  shell  theories  by  Naghdl  and 
Vangsarnplgoon  [89],  and  Llbrescu  and  Schmidt  [90].  In  all  these  works, 
no  attempt  was  made  to  obtain  solutions  of  the  theory. 

From  the  review  of  the  literature  on  various  plate  and  shell 
theories,  one  can  make  the  following  observations: 


1.  Analytical  solutions  of  refined  theories  of  anisotropic 
composite  laminates,  other  than  those  with  simply  supported 
edges,  are  not  available.  For  example,  the  Levy  type  solutions 
are  not  developed  for  laminated  composite  plates.  This  might 
be  possibly  due  to  the  difficulty  In  solving  the  ordinary 
differential  equations  resulting  from  the  Levy  procedure. 

2.  The  most  commonly  used  refined  theory  Is  the  first-order  shear 
deformation  theory,  commonly  referred  to  as  the  Mlndlln  plate 
theory.  It  Is  well-known  that  the  first-order  theory  (1) 
requires  the  introduction  of  shear  correction  factors  (to 
correct  for  the  constant  state  of  transverse  shear  stresses 
through  thickness),  (11)  disregards  the  effect  of  transverse 
and  normal  stress,  and  (111)  does  not  satisfy  the  stress-free 
boundary  conditions  on  bounding  surfaces.  It  Is  also  known 
that  the  determination  of  the  shear  correction  factors  for 
anisotropic  composite  plates  Is  not  established,  and  their 
values  depend  on  the  geometry  and  lamination  scheme  of  the 
laminate. 

3.  Use  of  refined  theories  with  moderate  and  large  rotations 
(situations  quite  commonly  encountered  in  the  analysis  of 
helicopter  blades,  turbine  blades,  etc.)  Is  not  Investigated, 
especially  for  anisotropic  composite  laminates. 

During  the  present  research  analytical  solutions  based  on  the 
space-time  concept  were  obtained  and  the  third  shear  deformation  with 
moderate  rotation  terms  was  developed.  A  first-ply  and  post-first-ply 
analysis  capability  based  on  the  first-order  shear  deformation  theory 
was  also  developed.  A  brief  summary  of  the  first  two  results  is 
presented  In  the  next  section. 


2.  BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


A  Third-Order  Theory 

Consider  a  laminated  plate  composed  of  N  orthotropic  layers, 
symmetrically  located  with  respect  to  the  midplane  of  the  laminate.  The 
governing  equations  of  the  refined  theory  are  based  on  the  following 
displacement  field  1 37-39 J ; 


“1  'u*zl*x-3<K>2<*x  +  ij>' 

u2  ■  v  .  2[*y  -  $  ♦  |S)| 


u3  =  w  (1) 

where  (u^.UgtUj)  are  the  displacements  along  the  x,  y  and  2  coordinates 

respectively,  (u,v,w)  are  the  corresponding  displacements  of  a  point  on 

the  midplane  of  the  laminate,  and  ip  and  *  are  the  rotations  of  a 

x  y 

transverse  normal  about  the  y-  and  x-axes,  respectively. 

The  cubic  variation  of  Uj  and  Uj>  through  laminate  thickness 
introduces  higher-order  resultants. 


Px  =  J  o,z*dz  (1  *  1,2,6) 
h  1 
'  2 


(R1,R2)  3  J  z^(a5,(J4)dz 


and  laminate  stiffnesses, 


<FirHij>  ■  S  „  qu<z  *z  >dz  <'-J 


-  1.2.6) 


■  v-.n  •  i,;"!*;*! 


‘  !  » Mi*. 


(°i j*F1  j)  * !  h  Oij<*  •*  ><**  -  4.5) 


For  symmetrical  cross-ply  laminated  plates,  the  following  stiffness 
coefficients  vanish  (39): 

^1J  *  ^1j  *  ^  ^or  *  1,2, 4, 5, 6 

*16  ’  *26  ‘  °16  *  °26  *  F16  '  F26  *  H16  ‘  H26  ’  0 
*45  "  °45  '  F45  *  0 

This  Implies  that  the  effect  of  coupling  between  stretching  and  bending 
vanishes.  For  such  laminates  the  governing  equations  are  given  by  (see 
137.38J): 

3  3  3  3 

21  »  *„  A  3 K  A  3  K  A  »  ♦„ 

— ?  *F11  — 3”  +  Hll^ - — 3  +  — 4)  +  F12  — 2*“  +  M12^ - ?)(' — 2*~ 

3hz  11  axJ  11  3hz  axJ  ax4  Li  axzay  u  3hz  axzay 


+ -4^) ♦ 

ax  ay‘ 


3  *x  4  3  *x  A  3 

— 2 —  +  H12^ - — 2 —  +  — 2 — +  F22  — 3 

ay  ax  3hc  ay  ax  axcayc  cc  ayJ 


3  3 

aJ*  a  v 


♦  H22^"  — 3^  +  +  2Fgg(  /  ♦  — JT”)  ♦  2H66^"  — T~ 

zz  3hz  ay-5  ay*  60  axzay  ayzax  3hz  ayzax 

3^*u  ?a4  »  .2  a*  .  a*  2 

+  zfii  *  IS7)I  ”  7  l°55(77 *  ^  * Fss(~ 

+  °44<“?  +  +  F44(-  +  "?)1  +  (A55(TF  +  “?> 

+  l?|)l+q,0  (2a) 

ay 

J11  +  °12  i*£  +  Fll(-  +  t?}  +  F«('  i?)(^ 
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J  ♦  *  *w  4  •  *- 

>  *  °M(-T  *  iliy'  *  F«*('  Z!)(TT  * 


•  K  » K 


3h  ay 


'  |A55(,x  *  m>  *  °55(’  ^)(,x  *  «>'  '  ^  ,fll 

*  Hn<-  r?)(rr  *  *  Fi2  TJii  *  Mi2<- 


r1  *  V- 


*x>  *  f55(‘  ?I)(* 


»s 

'  l*44(*y  *  lj>  *  °«4(-  ^><*y  *  SJ*1  ’  *  iju* 

2  2  2  2 
4  *  S  *  *v  2i K  »  *,  4  *  \ 

*  H66<-  ^  ^  *  FU  HW  *  H«(‘  17'W 

*  *1  )  *  Fjj  J  ♦  Hjjf-  "^j)( — J*  *  M[)l  *  Tf 

ax2  ay  22  ay2  22  3h2  ay2  ayJ  h2  **  ay 

*  V  *  %4<-  Jr«§  *  V  "°  <?c> 

Here  w  denotes  the  transverse  displacement,  #  and  *  are  the  rotations 

^  * 

of  the  normal  to  Midplane  about  the  y  and  x  axis,  respectively,  q  is  the 
distributed  transverse  load,  and  A^j,  D^j,  F^j,  Hjj  are  the  plate 


stiffnesses,  defined  by 


(3) 


h/2 


(01j*F1j‘H1j)  "  J  ?  Qi$)(*2.A*6)dz  (l.j  -  1,2,6) 


h/2 


(*0n  »2  4, 


(A1J*°1J*FU)  *  -  .Qii  (!•*  •*  )dz  (l.J  *  *.5) 


-h/2 


'1  j 


Here  q|*>  denote  Che  reduced  orthotropic  Moduli  of  the  k-th  lamina.  The 
boundary  conditions  of  the  refined  theory  are  of  the  fora:  specify 

w  or  Q„ 


>n 


♦.  or  "n 
V,  or  Nn, 


on 


where  r  is  the  boundary  of  the  nldplane  a  of  the  plate,  and 
Mn  -  Mjnjj  ♦  h2nj  ♦  2M6nxny 

N  «  (N,  -  N,)n  n  ♦  Mc(n2  -  n2) 
ns  '2  1'  x  y  6'  x  y' 

P  «  P.n2  ♦  P,n2  ♦  2Pen  n 
n  lx  2y  6xy 


(4) 


P  -  (P,  -  P.)n  n  ♦  P,(n2  -  n2) 
ns  '2  1'  x  y  6V  x  y ' 

A  »P„«  jP« 

Q  -  Q.n  ♦  Q,n  ♦  (— ^  ♦  ~) 

yn  yl  x  v2  y  3^2  v  as  an  ' 


(5) 


H.  -  p.  (1  -  1,2,6) 

1  1  3h*  1 

Ql  «  Q,  -  ^  «i  (1  -  1.2) 

h 

a  a  a  a  a  a 

an  nx  ax  +  ny  ay  *  as  nx  ay  ’  ny  ax 

The  stress  resultants  appearing  in  Eq.  (5)  can  be  expressed  In  teras  of 

the  generalized  displacements  (w.e  ,e  )  as: 

*  J 
11 


«l  ■  •**(%  *  *  °S5<-  *  5?) 


P1 

'  F11 

ax 

*F12 

ay 

*  Hu<- 

4  w  * 
3hz  ** 

'  ♦ 

A. 

+  H12^ 

-  -*?)( 
3hz 

ay 

ay 

a* 

*  a* 

.2 

A 

a* 

2 

f2 

’  F12 

X 

ax 

—1 

ay 

♦  "12<- 

4  w  a 

— V  'ax" 
3hz  ax 

l-  ♦ 

a  w» 

axe 

+  H22^ 

‘  “7^ 

3hz 

y 

ay 

♦M) 

ay 

3* 

A 

a* 

a* 

2 

*  F66( 

—1 

ax 

♦  — -] 
ay  1 

♦  H 

'«*<-  p><1T  * 

ax 

t  ♦  2 

U*L) 

axay' 

*1  *  °44<i?  *  «,) 


44' 


3hz 


'ay 


D55<l?  * 


*x>*F5  5<-^HS 


*x> 


(6) 


Analytical  Solutions 

The  present  study  deals  with  the  development  of  the  Levy  type  solu¬ 
tion  of  the  refined  shear  deformation  theory  of  Reddy  ( 12,131  for  sym¬ 
metric  rectangular  laminates  with  two  opposite  edges  simply  supported  and 
the  remaining  edges  subjected  to  a  combination  of  free,  simply  supported 
and  clamped  boundary  conditions.  The  state-space  concept  Is  used  to 


solve  the  ordinary  differential  equations  obtained  after  the  application 
of  the  Levy  solution  procedure. 

The  Levy  Method  can  be  used  to  solve  Eqs.  (2)  for  rectangular  plates 
for  which  two  opposite  edges  are  slnply  supported.  TSe  other  two  edges 
can  each  have  arbitrary  boundary  conditions.  Here  we  assu«e  that  the 
edges  parallel  to  the  y-axls  are  simply  supported,  and  the  origin  of  the 
coordinate  system  Is  taken  as  shown  In  Fig.  1.  The  slnply  supported 
boundary  conditions  can  be  satisfied  by  trigonometric  functions  In  x. 

The  resulting  ordinary  differential  equations  In  y  can  be  solved  using 
the  state- space  concept. 

Following  the  Levy  type  procedure,  we  assume  the  following 
representation  of  the  displacements  and  loading: 

» 

w(x,y)  -  £  H  (y)slnox 

m-1  " 


v  (x,y)  -  r  X  (y)cosax 
x  m»l  m 


*v(*.y)  ■  Z  Y_(y)s1nax 

y  m-1  " 


(7) 


q(x,y)  -  t  Q_(y)s1nax 

m-1  " 


where  a  -  and  Xm,  Y^  and  denote  amplitudes  of  w,  *x, 
and  q,  respectively.  Substituting  Eqs.  (7)  Into  Eqs.  (2),  we  obtain 

e.W"  ♦  e.W"  +  e,H  ♦  e.X"  ♦  ecX  ♦  e*Y'"  ♦  e,Y'  «■  Q  -  0 

in  c  m  j  m  4  n  j  it  o  h  f  n  m 


v;  *  v.  *  eiox;  *  'll-.  ♦  'l/;  ■ 0 


*i3w-'  •  *  *uw;  *  'is*;  *  'is*;  *  * 


(8) 


Where  primes  on  the  variables  Indicate  differentiation  with  respect  to  y, 

13 


and 


»,  ■  - 


*2  '  ?<^>V<H12  *  *W  *  *44  -  °44  *  $>’F44 


*3  *  -,|(^>I*Sl  *  (^)?F5S  "  jJ?  °55  *  *55' 


*  17  N  *  17  "»  '  "«  *  5?  M«' 

“3  3?  <F“  ‘  31?  H‘l)  *  “'7  ““  '  (i?)2f«  ‘  *55' 


(f?2  '  “Tf  H22^ 
3hZ  ZZ  3hZ  ZZ 


7  °*4  *  <7>?F44  *  *44 


“2  M12  '  *66  *  (H12  *  V  - 


*4  •  ‘9  '  *5 

°«  '  17  F«  *  (I7,2h« 


•?|-»U  *  17  Fll  '  (^?)?Hlll  %l  °55  •  (7)2fS5  •  *55 


*|012  *  °S6  '  <f12  *  F66>  *  (y)2(H12  *  "m11 


*6  *  *14  ’  *7  •  *15  "  "  *12 


°22  '  Z1  F22  "  (^)Hi 


!17  '  ‘,?l'066  *  11?  FW  •  *  7  °44  -  (lJ)Zf44  '  *44 
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Equations  (8)  can  b«  written  as: 


W'" '  -  c.W-  ♦  c,W_  ♦  cJ(.  ♦  c*Y'  +  c  Q_ 

■  in  Z  m  3  m  4m  o^i 


where 


*:  •  's-i  *  v.  *  v.  *  v; 

Yi  ■  Vi' '  *  c10“i  *  cllXi  *  C12V. 


u_  t  *4*6*17  ,V?.e  )/(t  +  V, 
*10  *10*16  *ie  2  1  ei< 


16 


16 


(10) 


c-  • 


(r-5  *  r-f^  •  e3)/(ei  *  r-> 
*10  *10*16  J  1  *16 


*11*4  +  !nVli  .  e  )/(e  +  !§_) 


10 


*10*16 


16 


,  ’6e17  .  *4*12  e6e12  r  Ulm  ,  e6  v 
(“5 —  ♦  -5 —  ♦  JT~ 2 ®7>  /  (®i  +  * — ) 

"i*  Cm  •  1  Cjg 


c:  *  - 


!16  *10  *10e16 
*16 

*1*16  *  *6 


C5  *  ■*4/*10*  c6  *  '*5/e10*  C7  *  -ell/e10*  c8  *  "e12/e10 

C9  *  *6/e16*  c10  *  *7/e16»  C11  *  e12/e16*  C12  *  'e17/e16  (11) 

The  linear  system  of  ordinary  differential  equations  (10)  with 
constant  coefficients  can  be  reduced  to  a  single  matrix  differential 
equation  using  the  state-space  concept  (see  (91]) 

x '  *  Ax  +  b  (12) 

This  can  be  done  by  Introducing  the  variables 

*i  ■  v  **  •  v  «3  •  *«  •  wi' 1 


*«•  x6  ’  *i’  *7  *  ¥»*  x 


-  Y ' 

8  M 


(13) 


where 
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01  0000  0  0  “ 

0  0  1  0  0  0  0  0 

0  0  0  1  0  0  0  0 

C2  0  CX  0  C3  0  0  C4 

A  ■ 

0  0  0  0  0  1  0  0 

c6  0  c5  0  C7  0  0  c8 

0  0  0  0  0  0  0  1 

_0  c10  0  Cg  0  cn  c12  0  _ 

The  solution  of  Eq.  (12)  Is  given  by 

«  ■  «AyK  ♦ 

where  K  Is  a  constant  vector  to  be  determined  from  the  boundary 
conditions,  eAy  denotes  the  product. 


eAy  -  (cl 


ex2*  0 


v  r  i’ 


(c)  Is  the  matrix  of  distinct  eigenvectors,  x^i  *  1,2, 3, ...8)  are  the 
eigenvalues  associated  with  matrix  A,  and  (cl-1  Is  the  inverse  of  the 
eigenvectors  matrix  [ c ] . 

The  following  boundary  conditions  are  used  on  the  remaining  two 
edges  (i.e.,  the  edges  parallel  to  the  x-axis)  at  y  =  * 

simply  supported:  w  =  ^  =  =  0 

clamped:  w  =  —  =  iti  =  =  0 

- E —  ay  x  y 


free: 


P2  =  M2  *  0 


- P-  =  0 

6  3h2  6 

<»*- "*  ♦£<**♦  **> 


(17) 


Numerical  results  are  presented  for  orthotropic  and  symmetric  cross- 
ply  (0°/90o/0°)  plates  subjected  to  three  types  of  loads:  uniformly 
distributed  load  (qQ),  triangular  distributed  load  (2q0)  and  concentrated 
load  P  as  shown  in  Fig.  2.  The  following  sets  of  material  properties  are 
used  in  the  calculations: 


Material  I: 


*  20.83  x  106  psi 
G^2  =  6.10  x  10®  psi 
G23  =  6.19  x  10®  psi 


,  E2  =  10.94  x  10®  psi 
,  G^  =  3.7  x  10®  psi 
,  =  0.44 


(18) 


17 


Material  II: 


Et  -  19.2  x  106  psl  ,  E2  •  1.56  x  106  psl 

g12  *  G13  "  0,82  x  1q6  psi  •  g23  *  °-523  x  1q6  ps1  (19) 

v^2  *  0*24 

The  following  notation  has  been  used  throughout  the  figures: 

SS  -  simply  supported  at  y  *  -b/2  and  at  y  »  b/2. 

CC  -  clamped  at  y  *  -b/2  and  at  y  *  b/2. 

FF  -  free  at  y  *  -b/2  and  at  y  ■  b/2. 

SC  -  simply  supported  at  y  *  -b/2  and  clamped  at  y  *  b/2. 

SF  -  simply  supported  at  y  *  -b/2  and  free  at  y  =  b/2. 

CF  -  clamped  at  y  *  -b/2  and  free  at  y  =  b/2. 

UN  -  uniformly  distributed  load. 

TR  -  triangular  distributed  load. 

PL  -  point  load  at  the  center  of  the  plate  (20) 

To  show  the  effect  of  transverse  shear  strains  on  the  deflections 

plots  of  nondl mens Iona 11 zed  center  deflection,  w  =  103w(a/2,0)h3E2/(qoa4), 
versus  side  to  thickness  ratio  of  various  plates  are  presented  In  Figs. 
3-5.  The  shear  deformation  effect  Is  more  significant  In  cross-ply 
plates  than  In  orthotropic  plates.  Also,  the  first  order  shear 
deformation  theory  (FSOT)  over  predicts  deflections  relative  to  the 
higher  order  theory  (HSDT). 

Figures  6  and  7  contain  plots  of  the  transverse  stresses  o^3 
through  laminate  thickness  for  various  boundary  conditions.  The 
stresses  were  computed  using  lamina  constitutive  relations.  The 
transverse  shear  stresses  are  constant  and  parabolic,  through  thickness 
of  each  lamina,  respectively,  for  the  first-  and  higher-order 
theories.  The  discontinuity  at  Interface  of  lamina  is  due  to  the 
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—  * -t„i /v*V  V  T’3  *  f  *4*  V  V  t1  *■'  *  *  V  V  fV ■  **  "tT 

.  -  i ,  y.  y,v>%y/-,y.y  vVv 


Ficrure  3. 


ionalized  center  deflection 
Mondinensi  SSSC  plates  using  theories 

»^f0Sraer  sV.aar  a.fom.t.on  <*«>“ 
and  hi^at  (HSW  oni£otin  load). 

(Materia-  2,  a/o 
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Figure  4.  Nondimensionalized  center  deflection  versus  side-to 
thickness  ratio  of  SSFC  plates  (Material  2,  a/fc=4). 


(0°/90°/0°),  FSDT 


a/h 


Nondinensionalized  center  deflection  versus  side-to 
thickness  ratio  of  SSCC  plates  (Material  2,  a/b  =4) 


Figure  6.  Variation  of  the  transverse  shear  stress  through  the  thickness  of  orthotropic 
plates  subjected  to  uniform  load  (Material  2,  a/fc*»4,  h/a=0.14). 


24 


Figure  7.  Vatidtion  ot  the  tiausverse  shear  stress  through  the  thickness  of  cross-ply  (0/°0/0) 
laminates  under  uniform  load  (Material  2,  a/b=4,  h/a-0.14). 


60 


0 


j _ i _ i _ l 

5  10  15  20 

a/h 


Fiqure  4  .  Norsdiner.s  lona  1 1  zed  renter  stress  versus  s  l  1e-to-th  lcxness 
ratio  for  SSFF  lanmates  under  umforn  load  (a  b*4. 

Mater l a  1  2 )  . 


mismatch  of  the  material  properties.  When  the  stresses  (°x*<Jy»oxy) 
obtained  from  the  constitutive  equations  are  substituted  Into  the 
equilibrium  equations  of  elasticity  and  Integrated  to  determine  the 
transverse  shear  stresses,  the  resulting  functions  will  be  continuous 
through  the  thickness. 

Plots  of  the  nondlmenslonal Ized  center  stress,  *  10?oll(f* 

h  2  2 

0,  j)h  /(q  a  ),  versus  side  to  thickness  ratio  for  simply  supported  and 
free-free  (SSFF)  plates  are  shown  In  Figs.  8  and  9.  The  shear 
deformation  effect  Is  quite  significant  for  a/h  ratios  smaller  than  10. 


A  Moderate-Rotation  Theory 

The  theory  Is  a  general Izatlon  of  the  classical  plate  theory,  the 
first-order  shear  deformation  plate  theory,  and  the  third-order  shear 
deformation  theories  of  Reddy  (37-39].  The  theory  Is  based  on  an 
assumed  displacement  field  and  orders  of  magnitudes  of  linear  strains 
and  rotations.  The  associated  strain-displacement  equations  are 
presented  and  the  equations  of  motion  are  derived  using  the  principle  of 
virtual  work.  Specialization  of  the  equations  of  motion  for  various 
existing  theories  is  presented. 

Points  of  a  three  dimensional  continuum  V  are  denoted  by  their 
orthogonal  curvilinear  coordinates  a  =  (xl,x2,x3).  Covariant  and 
contravariant  base  vectors  at  points  of  the  continuum  are  denoted 
by  ^  and  g1,  respectively.  Latin  Indices  are  assumed  to  have  values  1, 
2,  3,  and  the  Greek  indices  have  values  1,  2.  The  laminated  plate 
continuum  In  the  undeformed  configuration  Is  defined  by  the  Cartesian 
product  of  points  In  the  midplane  a  and  the  normal  [-  h/2,  h/2): 


where  h  denotes  the  constent  thickness  of  the  laminate.  Let  x°  denote 
the  curvilinear  Inplene  coordinates  and  x3  be  the  normal  to  a.  The 
■etrlc  tensor  components  of  q  are  denoted  by 

afl  a  a  33  • 

■  9.  •  s.  •  «  -  3  3  .  9  *  *33  •  1 

s.  ’  ^  •  9.  •  9*  ■  *!  •  i3  *  ;  <*» 

where  r  Is  the  position  vector  of  a  particle  (x°,x3)  at  time  t,  «fl  Is 

-  Ck 

the  Kronecker  delta,  and  n  Is  the  unit  normal  to  the  boundary  of  o. 

The  displacement  vector  of  a  point  In  the  plate  at  time  t  Is  of  the 

form 


u  *  uaga  ♦  u^n  »  u^g®  ♦  u^n 


where  the  Elnesteln  summation  convention  on  repeated  subscripts  Is 
assuewd.  The  covariant  components  of  the  Green- lagrartge  strain  tensor 
are  given  by 


■  ?  (utu  *  uj 1 1  *  u.nu"IJ) 


where  a  vertical  line  denotes  covariant  differentiation.  The  strain 
components  e i j  can  be  expressed  In  terms  of  the  linearized  strains  e^j 


and  rotations  as 


£1j  *  *1 j  +  1  em1ej  *  \  (em1wJ  +  emj“™)  +  ?  “ml“j  (24) 


m.  1  m 


where 


eij  "  \  (ui|j  +  UJM)  •  wij  m  \  (ui|j  ’  UJM) 

Following  [901,  we  now  assume  that  the  strains  and 
rotations  u. .  are  of  the  following  magnitude: 


tlj  -  0(a2)  .  u.aa  »  0(e2)  ,  «a3  -  0(e)  .  e  «  1  (26) 

Equation  (26)  Implies  that  the  strains  and  the  rotations  about  the 
noroal  to  the  el  dp  lane  are  small,  and  that  the  rotations  of  a  normal  to 
the  midplane  are  moderate.  Such  assumptions  are  justified  In  view  of 
the  large  Inplane  rigidity  and  transverse  flexibility  of  composite 
laminates. 

A 

Neglecting  terms  of  order  (e  )  and  higher  In  the  strain 
displacement  equations  (24),  we  obtain 

1  a  *  *  a  *  5  (*■>  ¥  e..  (I*3)  +  5 

afl  afl  2  3a  fl  36  a  2  3a  B 


ea3  +  \  (eXaW3  +  *33Wa)  +  \  WXa“3 


c33  *  e33  *  +  7  wx3w3  W 

where  the  underlined  terms  are  of  order  (e3). 

The  present  theory  Is  based  on  the  following  assumed  variation  of 
the  displacement  components  across  the  plate  thickness: 

u  (x6,x3,t)  ■  u°(xa,t)  -  x3u?,  +  f(x3)u1(xB,t) 

a  a  J  |  a  a 

u3(x8,x3.t)  -  u$(xfl.t)  ♦  u°(xB,t)  (28) 

where  f  Is  a  specified  function  of  the  thickness  coordinate  x^.  Note 
that  the  transverse  deflection  Is  assumed  to  be  Independent  of  x^  and 
consists  of  two  parts,  one  due  to  bending  and  the  other  due  to 
transverse  shear.  The  particular  form  of  displacement  field  Is  assumed 
In  order  to  Include  the  displacement  fields  of  the  classical  plate 
theory  (set  u^  ■  0  and  u*  *  0),  the  first-order  shear  deformation  theory 
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[set  Uj  *  0  and  f(x3)  »  x3l,  and  the  third-order  shear  deformation  theory 
of  Reddy  (391  (set  u^  -  0  and  f(x3)  «  x3[l  -  |  (jj-)2)l,  among  others. 


For  the  displacement  field  In  Eq.  (28),  the  strains  for  the 
moderate  rotation  theory  become  (consistent  with  the  assumptions  In  Eq. 
(26)1, 


e  *  e°  ♦  x3e^  ♦  ft 
a8  a8  aB  afl 


Ea3  *  ca3  *  qt°*3  *  *3"°3  +  9x3<l3  *  fS3  +  fg'a3 


o  -o  2*0 
c33  “  c33  ♦  g‘ 33  +  g  c33 


(29) 


where  g  »  df/dx3,  and 
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(30) 


The  dynamic  version  of  the  principle  of  virtual  displacements  Is 
used  to  derive  variational ly  consistent  equations  of  motion  associated 
with  the  displacement  field  In  Eq.  (28).  The  principle  can  be  stated. 
In  the  absence  of  body  forces  and  prescribed  tractions,  as 

0  -  J  If  <o1J5c.,)dV  +  /  q«u,dA  -  J  0(u,«u.)dVldt  (31) 
0  V  1J  Q  V  1  1 

where  a ^  denote  the  contravarlant  components  of  the  symmetric  stress 

tensor,  q  ■  q(xa)  Is  the  distributed  transverse  force  per  unit  area. 


and  0  Is  the  density  of  the  material  of  the  plate.  The  superposed  dot 
denotes  the  time  derivative,  u  =  au/at.  We  Introduce  the  couples  and 


Inertias, 


h 

(Nafl,  Hafl,  Pa8)  *  J  oa6(l,  x3,  f)dx3 

h 

‘  7 

h 

(Q°»  Q°,  Sa,  Sa)  =  J  aa3(l,  g,  x3,  x3g,  f,  fg)dx3 
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‘  2 


h 

(N3.  N3.  N3)  *  ;  o33(l,  g,  g2)dx3  (32) 
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*  2 

h  h  h 

2  2  2 

IQ  =  /  odx3  ,  Ij  =*  J  px3dx3  ,  l[  *  J  h  pfdx3 
h  h  -  -x 


h  h  h 

I-  =  J  p(x3)2dx3  ,  L  =  J  px3fdx3  ,  ij  -  J  Pf2dx3  (33) 
c  h  c  h  ^  h 

'2  '2  '2 

The  equations  of  motion  of  the  theory  are  obtained  by  substituting  Eq. 

(30)  for  the  strains  In  terms  of  the  displacements  (u°,u?,u?,u*) 

a  J  j  a 

Into  Eq.  (31),  Integrating  by  parts  to  transfer  differentiation  from  the 
displacements  to  the  stress  resultants  and  couples,  collecting  the 
coefficients  of  the  various  virtual  displacements,  and  invoking  the 
fundamental  lemma  of  the  calculus  of  variations.  We  obtain  the 
following  six  equations: 
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N°8lfl  -  «Sk>la  +  W^l,  *  ‘ou°  -  ‘ft|.  +  ‘ft 

ft=  H°,|a.MN"(u°|,^°|t)l|a-(V;|a)|. 

-  <R“i«ulia>ia  *  («“u*)ia«  -  (s°“iu>i.  +  (H3u;,  )i 


-  (nV)|  =  q  +  1  (u°  +  u°)  +  I.u0,  +  Lu1,  +  I9u° 
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*“r  <N°fl<u5|s  *  °?|s>Ha  -  Q“la  =  q  *  i0(u|  *  G§) 

ft=  pa9is  + 

-  <S8u3U>U  -  *  ^3la  ■  ‘ft  -  ‘ft|„  +  ‘ft  <34> 

where  the  underlined  terms  are  entirely  due  to  the  Inclusion  of  moderate 
rotations  (i.e.,  over  and  above  the  von  Karman  nonlinear  terms). 

Equations  (34)  can  be  specialized  to  the  three  different  theories 
discussed  earlier.  The  equations  are  summarized  below: 

(1)  Classical  Plate  Theory  (u?  =0,  u*  =  0) 


NaSl  -  (Q0u?i  )|  =  I  u° 
6  v  31a7  1 B  0  a 


M"“La  *  <N“ft|  A  -  «ftlA  -  ^IftlA  +  (“ftlJI 


=  q  +  I  u?  +  I,u0.  -  I-u°. 

0  3  1  a|a  2  3  |aa 

(ii)  First-Order  Shear  Deformation  Plate  Theon 
(u°  =  0,  f  =  x3) 


*  i  *i 
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M°aU  -  q8(5.a +  -  !!ft  *  Aft  -  ‘ft  ♦  ‘ft  <36> 
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(Ill)  Third-Order  Shear  Deformation  Plate  Theory 
(u°  -  0.  f  .  x3|l  -  |  (x3|h)2]) 
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Po8l  -  Q8(«  +  u°.  )  +  SB|  ul  -  NV  *  i;U°  +  I-u1  (37) 
'8  v  a8  a|  8  1 8  a  _ a  la  Z  a  v  ' 

Note  that  several  other  theories  can  be  obtained  from  Eq.  (35)  as 
special  cases.  For  example,  the  refined  theory  of  Kromm  (22,23),  along 
with  that  of  Basset  [5],  can  be  obtained  by  setting  *  0. 


c&i  i.l  ki3.,l 


rf::o 


rf.“.l 
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